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^ ■ Abstract 

' Nucleosynthesis restrictions on mixing of active neutrinos with possible sterile ones 

. are obtained with the account of experimentally determined mixing between all active 

neutrinos. The earlier derived bounds, valid in the absence of active-active mixing, are 

reanalyzed and significant difference is found in the resonance case. The results are 

obtained both analytically and numerically by solution of complete system of integro- 

QQ . differential kinetic equations. A good agreement between analytical and numerical 

<^ ' approaches is demonstrated. A role of possibly large cosmological lepton asymmetry is 

95 . discussed. 
O 
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Combined data from SNO [1], KamLAND [2], and earlier data on solar [3] and atmo- 
Q-i! spheric [4] neutrino deficit present strong indication that all these neutrino anomalies can 

^ ' be explained without introduction of a new sterile neutrino, for recent analysis see refs. [5]- 

[11]). According to ref. [8], a mixing of active, Va (a = e,fj,,T), and sterile, Ug, neutrinos 
is excluded at the level sinr; < 0.13. On the other hand, at a weaker level, mixing with 
^ ■ Vg is not only allowed but even desirable. According the ref. [12], some features of the 

^ . solar neutrino data are better described if there exists fourth sterile neutrino state with 

the mixing parameter sin^ 29 = 10~^ — 10"'^ and the mass difference Sm"^ = (10~^ — 10~^) 
eV^. Moreover, the LSND results [13] probably also demands mixing with additional one 
or several new (sterile) neutrinos or maybe, if confirmed, something even more drastic. 

Consideration of big bang nucleosynthesis (BBN) may noticeably strengthen the upper 
bound on mixing between active and sterile neutrinos and even exclude explanation of 
LSND anomaly by admixture of Ug. Such limits can also restrict explanations of solar and 
atmospheric neutrino features by an admixture of Ug. BBN bounds on mixing between Ug 
and z/(j were derived long time ago [14]- [19] under assumption that only one active neutrino 
is mixed with a sterile partner and mixing between active neutrinos was disregarded (for 
more references and discussion see the review [20]). A simplified consideration of 4 neutrino 
mixing has been done recently in ref. [21] for rather large mixing angle, sin^ 29ag > 0.01. 

It is practically established now that all active neutrinos, u^, Vfi, and v^-, are strongly 
mixed with the parameters given by Large Mixing Angle (LMA) solution to solar neutrino 
deficit [11]: 

5m'^^o\ = 7.3 • 10"^ eV^, tan^ O^^i = 0.4 (1) 
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and by atmospheric neutrino data [4]: 

Sm\tmo = 2.5 • 10-3 eV^ tan^ ^^tmo ~ 1 (2) 

Existence of fast transitions between z/g, i^^, and v-r noticeably changes BBN bound on 
mixing with sterile neutrinos especially for small values of mass difference. The reason for 
this change is the following. There are three possible effects on BBN created by mixing 
between active and sterile neutrinos. First is the production of additional neutrino species 
in the primeval plasma. The second effect is a depiction of the number density of electronic 
neutrinos which results in a higher neutron freezing temperature. Both these effects lead to 
a larger ncutron-to-proton ratio and to more abundant production of primordial deuterium 
and helium-4 (for the details see e.g. review [20]). If mixing between active neutrinos is 
absent then the second effect would manifest itself only in the case of (fg ~ ^'s)-mixing, if 
we neglect relatively weak depopulation of through the annihilation i/gVe — ^ i'ij.,tI'ij,,t- In 
the realistic case of experimentally measured mixing between all active neutrinos a deficit 
of Vfj_ or v-j- would be efficiently transformed into deficit of leading to a stronger BBN 
bound on active-sterile mixing. 

The third effect is a generation of large lepton asymmetry due to oscillations between 
active and sterile species [22]. However, this effect takes place only for very weak mixing, 
much smaller than the experimental bound and is neglected throughout almost all the paper, 
except for sec. 7. 

The paper is organized as follows. In sec. 2 general kinetic equations for oscillating 
neutrinos in the primeval cosmological plasma are presented. In sec. 3 the approximations 
that we use for analytical and numerical calculations are discussed and the effects of neutrino 
on Big Band Nucleosynthesis (BBN) are described. In sec. 4 we reconsider the case of mixing 
between one active and one sterile neutrinos, assuming that other active neutrinos are not 
mixed. Our calculations agree with the earlier works in non-resonance case, while there is 
a noticeable difference in the resonance case. In. sec 5 we consider for illustrative purposes 
algebraically much simpler situation when only two active neutrinos are mixed. The realistic 
case of all three active neutrinos mixed is presented in sec. 6. In all previous sections the 
calculations are done for negligibly small cosmological lepton asymmetry. The impact of 
the latter is discussed in sec. 7. In conclusion we summarize and discuss the derived bounds. 



2 Kinetic equations 

We assume that in addition to three known active neutrinos there exists a fourth neutrino 
flavor which does not have any interactions except for mixing with active neutrinos. The 
transformation between flavor eigenstates Ua and mass eigenstates Uj is described by the 
orthogonal matrix 

= Uaj Vj (3) 

where a = e,fi,T,,s and j = 1,2,3,4. We disregard here CP-violation effects, so matrix U 
is assumed to be real. In the limit of small mixing i/g ~ i^i, ~ 1^2, z^t ~ ^^3) and z/^ ~ 1/4. 
We know, however, that the real mixing between active neutrinos is large and they cannot 
be taken clS 3f dominant single mass cigcnstate. 

Kinetic equations for the density matrix of oscillating neutrinos have the usual form [23] : 

ip=[H^'\p]-i{H^^\p} (4) 

where the flrst commutator term includes vacuum Haniiltonian and effective potential of 
neutrinos in medium calculated in the flrst order in Fermi coupling constant, Gp. The 
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effective potential is proportional to the deviation of neutrino refraction index from unity 
and is calculated in ref. [24]. The second anti-commutator term includes imaginary part of 
the Hamiltonian calculated in the second order in Gp ^. This term describes breaking of 
coherence induced by neutrino scattering and annihilation as well as neutrino production 
by collisions in primeval plasma. 

Written in components kinetic equations have the form: 

ipss = ^'H^sKpas- Psa) (5) 

a 

ipaa = 'H^ah (Pba " Pab) + 'X-^a) (psa - Pas) - Hcollip) (6) 

b 

ipab = T^^Jpsb - 'H^sbPas + {^acPcb - '^a Pac) - ijabPab (7) 

c 

iPsa = 'Hi^Jpsa - 'Hi^Jpss + Yl i^^sb Pba " '^bJ Psb) " HasPsa (8) 

b 

where a, b, c label active neutrino species. Matrix elements of the first order Hamiltonian 
~ (below we will omit the upper index 1) are given by: 

Hap = Fa/3 + E irny2E)UajUpj (9) 

3=1-4 

where E is the neutrino energy and a and f3 run over all neutrino flavors, active and sterile. 
The effective potential in the medium, V^^, vanishes if any a or P are equal to s. For active 
neutrinos diagonal components of the potential are given by [24]: 

F„„ = ±Ci7?WGFr=^ - C^^^^ (10) 

a 

where T is the plasma temperature, Gp is the Fermi coupling constant, a = 1/137 is the 
fine structure constant, and the signs "it" refer to neutrinos and anti-neutrinos respectively. 
The first term arises due to charge asymmetry in the primeval plasma, while the second 
one comes from non-locality of weak interactions associated with the exchange of or Z 
bosons. According to ref. [24] the coefficients Cj are: Ci « 0.95, ~ 0.61 and Gl^'^ « 0.17 
(for T < rrifj). These values are true in the limit of thermal equilibrium, otherwise these 
coefficients are some integrals from the distribution functions over momenta. The charge 
asymmetry of plasma is described by the coefficients 77^") which are equal to 

n^"^ = ^Vi^e + Vui. + Vi^r +Ve- r?„/2 (for Ue) , (11) 
V^'"'^ = 2r/,^ + r?,, + - Vn/2 (for i/^) , (12) 

and Tj^'^^ for v^- is obtained from eq. (12) by the interchange fi t. The individual charge 
asymmetries, ijx, are defined as the ratio of the difference between particle-antiparticle 
number densities to the number density of photons with the account of the 11/4-factor 
emerging from e'''e^-annihilation: 

Tlx = (Nx - Nji) /N^ (13) 

If z/z/-interactions are essential then off-diagonal components Vab are non-vanishing [25] 
(see also recent discussion in refs. [26]- [28]). These components are proportional to the 

^Here and in the following we will use a natural system of unit in which % = c = 1. 
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integrals over neutrino momenta from ofF-diagonal components of the density matrix. In 
the case under consideration these non-diagonal components of the effective potential are 

sub-dominant (see sec. 5) and will be neglected. 

The coherence breaking terms in off-diagonal components of the density matrix are 
given by 7^/3 = (7^ + ^fp)/2^ where 7s = and 7a is the total reaction rate, including elastic 
scattering and annihilation of Ua- 

la = 9a GpT p . (14) 

Here p is the neutrino momentum and the coefficients ga, according to ref. [19] arc g^^ ~ 4 
and gu^,^j.r — 2-9. Slightly smaller results g^^ ~ 3.56 and gv^,ij,.r — 2.5 were obtained in 
ref. [29] due to account of Fermi statistics. 

The coherence breaking terms in kinetic equations for diagonal components of density 
matrix, paa, are given by 

holl = 2^ X! / '^'^(^2, 1^3, k)\Ael\^ [Paa{Pl)fl{P2) - Paaipz) fliP'i)] - 

2^ H / "^"^(^S, I'a, ^4)|^onn| [Paa{pi)Paa{p2) - flipzlflipi)] (15) 

where Agi^ann are the amplitudes of the corresponding reactions and the sum is taken over 
all possible channels of elastic scattering (first term) and annihilation (second term). The 
phase space element is given by: 

dr{k,u^M) = {2n) S {p, +P2-PS-P.) (16) 

It is noteworthy that in the absence of annihilation Ua^a ^ e^e~ the total number 
density of all neutrinos is conserved in comoving volume: 



a\t) j 



3ss + ^Paa] = COnst (17) 



(27r) 



where a{t) is the cosmological scale factor. At the temperatures below freezing of the 
reactions UaVa ^ e+e" (T/^ « 3.2 MeV and T/^ w 5.3 MeV [20, 30], see also sec. 3.2) 
sterile neutrino states are produced at the expense of active ones and the effective number 
of additional neutrino species does not change. However the effect of the oscillations on 
BBN would be noticeable even in this case because of the deficit of induced either by 
direct mixing — Vs or by mixing of z^^^,- — Ug and fast redistribution of active neutrino 
species due to large mixing between them. 

Concluding this section, we introduce convenient variables in terms of which we will 
solve kinetic equations: 

X = ma{t) = m/T, y = Ea{t) = E/T (18) 

where m is an arbitrary normalization mass which for convenience is taken as 1 MeV and 
a{t) is the cosmological scale factor. We assume that the temperature evolves as T = l/a{t). 
Though it is not always true, the precision of this assumption is sufficiently good for our 
pTirposcs. In terms of these variables the l.h.s. of kinetic equations in Priedman-Robertson- 
Walker background can be rewritten as 

f = {dt- Hpdp) f = Hxd^f (19) 
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where H = a/a is the Hubble parameter. 

Expressed through x and y, essential quantities take the following forms: 

H = h/x'^, with h = 4.46 • 10^22 (^^/io.75)^/^ (20) 

where g^^ is the number of species in the primeval plasma; g^^ = 10.75 for T < m^] we 
have neglected here presumably small (< 1) contribution from Vg. The energy difference of 
neutrinos is given by 

5E^^-^ = h- 10-13 (Sm'/eW^) {x/y) = dm{x/y) (21) 

Below we will always measure neutrino mass difference in eV^ if not stated otherwise. The 
"non-local" contribution to neutrino potential in matter (the second term in eq. (10)) is 
equal to 

Vi^:^ = C^''\y/x') (22) 

where C^^^ = 1.137 • lO'^o and ci^^ = c[^^ = 0.317 • lO^^o. Contribution to the potential 
from charge asymmetry of primeval plasma is always sub-dominant if r] has "normal value 
10"^ — 10~i°. The case of non-negligible charge asymmetry is considered in sec. 7. The 
coherence breaking terms determined by 7 have the same dependence on x and y as V^f^ but 
with much smaller coefficient. Still they should be retained because only they contribute 
into imaginary parts of the density matrix and create exponential damping. They are given 
by: 

lab = eab{y/x^) (23) 
where Cee = 1-5 • 10"^^ and e^a = ^tt = 1-1 • lO"^^. 



3 Preliminary considerations 

In this section we introduce some approximations and discuss some features of cosmological 
evolution of neutrinos which can be helpful for determination of numerical and analytical 
solutions of kinetic equations and for description of the effects of neutrino oscillations on 
BBN. 

3.1 Stationciry point approximation 

Brute force numerical solution of integro-differential kinetic equations (5-8) is rather difficult 
and for large mass differences suffers from numerical instabilities, especially in the resonance 
case. On the other hand, for small mass differences numerical calculations can be rather 
accurately performed. 

In order to avoid numerical problems, in this work we have used an approximate but 
rather accurate method discussed in ref. [30, 31]. If damping terms, given by lapPafi^ are 
sufficiently large, differential equations for off-diagonal components of density matrix can 
be formally solved in the stationary point approximation. It simply means that r.h.s. of 
these equations is assumed to be equal to zero and they reduce to linear algebraic equations 
(see eqs. (33,34) below). Such approximation works well when reaction rates are large, 
much higher than expansion rate. The latter is true for sufficiently large mass differences 
of oscillating neutrinos. It can be shown [14] (and it follows from the calculations presented 
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below) that in non-resonance case, the maximum production rate of sterile neutrinos takes 
place at the temperature 

r;;„^ = (10.8, 13.4) (3/y)i/3 (cos 2ef'\5m^ /eY^f'^ MeV (24) 

The first number above is for mixing of Vg with u^, while the second one is for mixing with 
or Vr- Thus for cos 29 5m?' > 10~^ eV^ the production of sterile neutrinos is efficient 
when T ^ H and stationary point approximation is sufficiently accurate. 

In this approximation one can express algebraically all off-diagonal components of (4 x 4) 
density matrix of oscillating neutrinos through 4 diagonal components, paa (a = e, /i, r) and 
pss- This is true if the off-diagonal components of effective potential of active neutrinos 
can be neglected (see sec. 5). Found in this way expressions for off-diagonal components 
of the density matrix can be inserted into eqs. (5-6) and we obtain a closed system of 
integro-differential equations for the diagonal components. The latter can be easily solved 
numerically as is done in the following sections. 

We remark that, while we have used quasi-stationary approximation for massive calcu- 
lations, we have checked the validity of this approximation by solving numerically the whole 
system of equations (i.e. without using the stationary point approximation) for several rep- 
resentative values of the parameters. This more exact approach takes much more computer 
time but we have not encountered any problem otherwise. These more accurate results are 
always in a very good agreement with those obtained with the approximate code. 



3.2 Kinetic equilibrium approximation 

Solution of kinetic equations is conceptually simpler in the case of a large mass difference 
between neutrinos. In this case production of predominantly takes place at the tem- 
peratures (24) which are much larger than the neutrino decoupling temperature (situation 
may be different for the resonance case, see below sec. 4). Thus at least for non-resonance 
transition the distributions of active neutrinos should be close to the equilibrium one. 

We will estimate neutrino decoupling temperatures in the Boltzmann approximation, 
i.e. we will assume /eg = exp{—E/T) = exp(— y). Under this assumption, if we take into 
account only direct reaction terms in the collision integral, the evolution of the distribution 
of non-oscillating neutrinos is governed by the equation: 

where the coefficient D is equal to 80(1 + g\ + 9%) if all reactions in which neutrino can 
participate are taken into account. Here ql = ±1/2 + sin^ 9w and gR = sin^^vi^, plus or 
minus in gi^ stand respectively for or i^h,t- The weak mixing angle 6w is experimentally 
determined as sin^ 9w = 0.23. Correspondingly the decoupling temperatures with respect 
to the total reaction rate are T,y^ = 1.34 MeV and r^,^,t/^ = 1.5 MeV. For this and the 
following simplified estimates the thermally average value (y) = 3 was taken. 

On the other hand, neutrino re-population, as one can see from eq. (6), is determined 
by the much weaker (inverse) annihilation rate. As a consequence, in some range of temper- 
atures, active neutrinos may be in kinetic but not in chemical equilibrium. In Boltzmann 
approximation this corresponds to: 

Paa{x,y) =na{x)eyip{-y) a = e,i2,T (26) 

In particular, if annihilation of active neutrinos into e'^e^-pairs is switched-off the total 
number density of active plus sterile neutrinos is conserved. The rate of annihilation I'v ^ 
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e+e" is given by eq. (25) with D = 16(g^ ^iid the decoupling temperatures would be 

T^^ = 3.2 MeV and T^^,^^ = 5.34 MeV. Below these temperatures the total number density 
of all four neutrino flavors in comoving volume could not change. So extra neutrino species 
are not created, though redistribution between Va and i^s is possible. 

In order to keep the computation time reasonably low, we assumed that, in the consid- 
ered parameter range, the active neutrino distribution functions are conveniently described 
by eq. (26). For we did not made any simplifying assumption and took it as an arbi- 
trary function of time and energy. Let us remark that, while we used kinetic equilibrium 
approximation for massive calculations, we have also solved numerically, for several values 
of the parameters, the whole system of integro-differential equations (i.e. assuming neither 
kinetic equilibrium nor stationary point approximations). The results obtained are always 
in a very good agreement with those obtained by using the approximate code. 



3.3 Effects on BBN 

One should also keep in mind that BBN is especially sensitive to the number density of 
electronic neutrinos because the latter directly participate in reactions of neutron-proton 
transformation. A depletion of and equally number densities ^ through oscillations 
would give rise to a larger neutron freezing temperature and to a larger abundances of pri- 
mordial deuterium and helium-4, thus enhancing the effect of the Vg contribution into the 
energy density of the universe. The described effect can be simply taken into account, if we 
assume that electron neutrinos are in kinetic equilibrium, that the Boltzmann statistics is 
valid, and neglect the electron mass. In this assumption, the neutron to proton intercon- 
vcrsion rate Tnp 

is proportional to F^p = (1 -|- ne)/2. This means that the neutron freezing 
temperature scales as: 



1 /2 

9* 



(l + ne)/2 



1/3 



(27) 



where 5* = 10.75 -|- {7/4:)ANy, and thus the global effect on BBN can be described as a 
variation of the effective number of neutrino species: 
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10.75 + (7/4) AAT, 



((l + ne)/2)2 



(28) 



If we consider possible deviations from kinetic equilibrium, the same formulas apply but 
we have to replace Ug by: 

, V _ /o°° dyy'^jy + Am/Tf Peeix,y) 

!^dyy^{y + Am/TYe-y ' ^^^^ 

where Am = 1.3 MeV is the neutron-proton mass difference and T is the temperature. 

If only one active and one sterile neutrinos were mixed, while mixing among the active 
neutrinos was absent the described effect would lead to a stronger restriction on the mixing 
angle in low dm? region in the case of (i'e-i^s)-™ixing. In the case of or v-j- mixing with 
Vg the effect is less pronounced because depletion of fg is now a two step process. First, 
the number density of or v-,- decreases due to transition into and after that some 
re-population of u^^t by Pg'^e ^ij,,t^ij,,t occurs at the expense of the number density of u^. 

^In the resonance case Ue and 9^ can be depleted by a different amount and a noticeable charge asymmetry 
between active neutrinos can be created. This is discussed in sec. 7 
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This process however is not very efficient because to be noticeable it should take place after 
decoupling of neutrinos from e+e" plasma but practically simultaneously neutrinos decouple 
from themselves. However if the active neutrino species are mixed then the depopulation 
of fe would be much stronger as we see in what follows. 



4 Mixed one active and one sterile neutrinos 

As a first step let us consider the case when mutual mixing of active neutrinos is absent 
and only one active neutrino is mixed with the sterile neutrino. Previous BBN constraints 
on sterile neutrino admixture have been derived under this assumption. We present below 
the results of our numerical calculations together with some useful analytic estimates both 
for the non-resonance and resonance cases. 



4.1 Non- resonance case 

The density matrix of the oscillating neutrinos is 2 x 2. It contains only four real functions: 
Paaj Pssj SLud Pas = Pas = R + H , which have to be determined from the solution of kinetic 
equations: 

Hxdxpaa = iT-tas {Pas - Psa) - holM) (30) 

HxdxPss = -iUas{Pas- Psa) (31) 

Hx dxPas = -i [{T-iaa - Ugs) - ilas] Pas + iUas {Paa - Pss) (32) 

The formal solution of equation (32) is: 



Pa 



f r ^as{Paa - Pss)l ■ f j i^aa " T^-ss - IJas) 

/ dxi YTr~\ / "^2 TTTl 

Jo {Hx)i J XI {Hx)2 



(33) 



where the indices sub-1 and sub-2 indicate that the corresponding expressions are taken at 
xi or X2- If 7as is large (more precisely we require 'jas/H ^ 1) then the integrals "sit" on 
the upper limit xi = X2 = x and in this way the stationary point result is obtained: 

Pas = -pr. ^^^T : (Paa - Pss) ■ (34) 

This result is valid both for resonance and non-resonance cases if the mass difference is 
sufficiently large, but the limits of applicability are somewhat different. As is mentioned 
above, for non-resonance case the temperature of Vs production (24) should be higher than 
the decoupling temperature with respect to total reaction rate. This condition is satisfied 
for if (5m^ > 3.7 • 10~^ eV^ and for v^^^t if bvr? > 2 • 10~^ eV^. The resonance case is 
considered in the following section. 

Substituting the result (34) into eq. (31) we find: 

Hx dxPss = 2 ^"C^''^ 2 (/^- - Pss) (35) 

which, in terms of the mass difference, dm?, and of the vacuum mixing angle ^ , 9, can be 
rewritten as: 

HxdxPss = -r [Paa - Pss) — ,.^.2 , ^^^^ 

4 {cos2e -Vaa SE) ^ii ^bE^ 



^In our notations, is positive if sterile neutrino is heavier than active neutrino, in the Umit Q 
Non resonance cases correspond to positive 5m^. 
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Prom this expression (neglecting the 7^-term in the denominator) follows that the rate of 
production of sterile neutrinos in the primeval plasma is equal to: 

7. = ^7atan2 2^„ (37) 

where 9m is the effective mixing angle in matter given by the last factor in eq. (36). This 
result is twice smaller than simple estimates presented in earlier papers [14]- [18]. 

In non-resonance case, when cos 29dE — Vaa / (i.e. 5rn? > 0) and one may neglect 
the 7^-term in the denominator, the number of additional neutrino species ANi, at BBN 
can be analytically calculated from eq. (36) in the limit of pss <C Paa and paa = feq- In this 
approximation we obtain: 

((5m^,^,/eV2) sin^ 29"^"^ = 3.16 • IQ-^AN^f (38) 
(Sm'^u^uJeV^) sin^ 20^'^^^ = 1.74 • lO'^iAN^f (39) 

If ANi, is not small, then pgg in eq. (36) should be included and a better approximation in 
the bounds above would be given by the factor ln^(l — AA'^^) instead of (AN^)"^. 

The described results are noticeably weaker than earlier obtained analytical bounds 
They well agree, instead, with our numerical calculations and with those of ref. [19]. Our 
numerical results are presented in figs. 1 {vi^ — Vg mixing) and 2 (z/g — Vs mixing). In the 
first panel, we show the effects of oscillations on the energy density of Uf.- Specifically, wc 
show the lines which correspond in the plane {sin^ 29 , dm^) to fixed, but different, values 
of the energy density of electron neutrinos (normalized to the equilibrium value) ^. In the 
second panel, we show the energy density of sterile neutrinos. In the third panel we show 
the total neutrino contribution to the energy density, each line corresponding to a certain 
value of the number of extra neutrino flavors AN^. If there was not an additional impact of 
Vf, on the neutron-proton transformations, these lines would present BBN bounds on mixing 
parameters. The total effect of oscillations on BBN, including the impact of on n — p 
transformation according to eq. (28), is demonstrated in the fourth panel in terms of the 
effective number of neutrino species AA^^^'"^. Finally, we show with red dotted lines the 
approximate analytical results described above. Numerical calculations were made by using 
stationary point and kinetic equilibrium (for active neutrinos) approximation. We checked 
by comparing with the complete calculations, that these approximations do not introduce 
significant errors. 

As one can see, the bounds (38,39) are reasonably accurate for large mass differences 
since, in this case, re-population of active neutrinos is very efficient and we can always 
assume paa = feq- On the contrary, when Sm'^ is small (say, Sm? < 10~^ eV^) significant 
deviations are expected. In this case, sterile neutrino production takes place at low tem- 
peratures, see eq. (24), when inverse annihilation processes become slow with respect to the 
expansion rate of the universe. Slow re-population of active neutrinos results in a decrease 
of AN^, since the production of Ug goes at the expense of diminishing the number density 

of l^a- 

In the case of Ve — mixing, the decrease in the number density of and Ue results in 
a decrease of the neutron-to-proton inter-conversion rates which leads to a higher temper- 
ature of neutron freezing and to larger neutron-to-proton ratio and, as a result, to higher 
abundances of D and "^He. This effect can be effectively described as an increase in 
which is much larger than the simultaneous decrease of AN^ produced by the decrease of 

''Numerical calculations have been mostly done in the assumption of kinetic equilibrium for active neu- 
trinos, as explained in sect. 3.3. In this assumption, the energy density of i>e (normalized to equilibrium 
value) is simply equal to the value of the overall normalization factor Ue defined in eq. (26). 
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fg number density. This means that, for small Sm^, the BBN bound on Ug — Ug oscillations 
is considerably stronger than that obtained from the simple estimate given by eq.(38). 

It is important to remark that in the case of — [v^- — Vg) mixing, one expects a 
depletion of too. The number density of (or u^^) goes down due to transformation of 
these particles into Vg. This opens window for Vf.Ve. l^^l,Tl^^l,T and to a decrease of number 
density of u^,- In this case, however, the effect is considerably smaller and it has a magnitude 
comparable to the the decrease of AA^,^ produced by the decrease of {vr) number density. 



4.2 Resonance case. Analytic estimates 

Situation in the resonance case (i.e. when 5w? < 0) is somewhat more complicated both 
from the point of view of analytical and numerical calculations. We show, in this section, 
that simple analytical estimates can be obtained even in this case. 

Let us return to kinetic equations (30,31). We assume that the collision integral Icoll{l) 
has a negligible role during resonance. In other words, we assume that during resonance 
neutrinos are neither effectively created nor scattered by other particles. This approximation 
is justified because the characteristic "time" of collisions/annihilation, 5xcoii ~ l/H^ is much 
longer than the inverse resonance width, Sxyes ^ l/^E, in the essential range of parameters. 
We maintain, however, the decoherence term in the equations for the non-diagonal terms 
of the density matrix. Under this assumptions, neutrino evolution during resonance is 
described by the equations: 

Hxd,{paa + Pss) = (40) 
Hx dx{paa - Pss) = 2i Has {Pas - Psa) (41) 

where pas = pla the non diagonal term p^^ is given by eq.(33). We can write formally: 

Pas{x,y) = ieM-n^,y)) /"(ia;il^^^^^%^^^exp(F(xi,y)) (42) 

Jo \Hx)\ 

where the function F{x, y) is defined by the condition: 

Hx d^F{x, y) = i (Haa - Ugg - i-fas) (43) 

From the above expression, we see that the function F{x, y) has stationary points Zres in 
the complex plane where: 

Uaa {Zres,y)- Ugs {Zres,y) - ilas (^res , y ) = (44) 

We can explicitly solve this equation, obtaining: 



1/6 / \l/6. 



( ^ _ iy'{\Ci\+^eas)\ ( rlQI 



\dm\cos{26) I \\dm\cos{29) 







exp (^i^^ (45) 



where n = 0,1,. ..,5, the parameters dm, and eas are defined in eqs. (21)-(23) and 
5a = \eas/Cy \ ~ 10~^. One can see that the usual resonance condition, Haaixres^v) — 
'Hss{xTes,y) = 0, gives the resonance point Xr 



^res- 



\\dm \ cos(20) / 
which is essentially equal to Xres — Re ((-2;res)o)- 



y. 1 (46) 
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We can take the integral (42) by using the saddle point approach: 



Pas{x,y) 



iV27r exp [-F{x, y) + F{Zres, V)] 0{x - Xres) 

-1/2 



'Has 






[{Hx)\ 


res 





exp \ «^ - ) {poAx - Pss)x 



(47) 



)o, and 



where the index [res] indicate that the various quantities arc evaluated at Zrcs = ( 
the phase ^ is given by = arg(5^F/(?,T'^)i.cs- Wc remark that the second derivative of the 
function -F(x, y) at resonance, can be expressed througii the simple formula: 



d''F{x,y) {n 



vac 
aa 



' ^ss Jres 



(Hx) 



(48) 



where TC^^^ is the vacuum hamiltonian, from which wc sec that -0 ~ tt/2. 

Equation (47) describes the approximate behavior of pas{x,y) as a function of neutrino 
parameters and of the difference (pao — Pss)res at the resonance. We can use this expression 
to integrate eqs. (40,41). We have: 



(Pao + Pss) 
MPaa - Pss) 



= const 

rx 

= 2i 

Jo 



dxo .rj \ \Pa 



PsajQ 



(49) 

(50) 



where IS.{paa — Pss) is the total variation of {paa — Pss) across resonance. If we integrate the 
previous equations using, one more time, the saddle point approach, we obtain: 



d^F 


-1 - 


dx'^ 


res - 



^{Paa - Pss) = -27r(paa - Pss)res 

This expression can be calculated explicitly, giving: 

_ TT \d^\sm\2e) 
l^KPon Pss) - g hcos{2d) 

TT \d^\9.\^'^('?.Q\ 

Re 



{'Has)! 



(Hx), 



exp [i— — iip ] + c.c. 



(51) 



Rel^l ipaa-Pss)r 



6 ^cos(2^) 



|d^|cos(20)) 



1/2 



{Paa - Pss) r 



(52) 



As a final result, we have thus obtained an algebraic relation between the value of the 
difference {paa — Pss) at the resonance, {paa ~ Pss)tcs, and its total variation A{paa — Pss)- 
If we assume that A{paa — Pss) is small, the value of (paa — Pss) at resonance is close 
to its initial value, i.e. {paa — Pss)res — {Paa)in- In this assumption, we immediately obtain 
from eqs. (52,49) an estimate for active (sterile) neutrino depletion (production) due to 
resonance. We have: 



P — 

^ aa — 



P = 
^ as — 



{pa 



^Pa 



{Paa) 



1 - K 



Aps 
{Paa) 



where 



^ TT |d„|sin2(20) fymcS\+ie, 
ri = — — : — Ke 



= K 



1/2 



12 /icos(26l) 



\dm\ cos (2^) 



Aa 



1/2 



sin^(2g) 

eV^ ) cos(2e)3/2 



(53) 
(54) 

(55) 
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The coefficient is equal to Ag = 4.4 • 10'^ for i/g — i^s mixing and A^^^ = 2.3 • 10^ for — Vg 
or Vr — Vs mixing. The previous equations represent our final results. A few comments are 
in order: 

1. These simple formulas, in the range of parameters in which they are valid (i.e. small 
K and, thus, small mixing angles), well agree with the results obtained by describing 
resonance effects by Landau-Zener approach [39]. We recall that, in this scheme, the 
active (sterile) neutrino survival (oscillation) probability Paa {Pas) is given by: 

Paa = Pc (56) 
Pas = 1-Pc (57) 

where: 



Pe = exp i-K) (58) 

and we assumed sin 2^ <^ 1. The agreement between our simple formulas, which 
take into account decoherence effects, with the Landau Zener results (which do not 
take them into account), in the range of parameter in which our formulas are valid, 
makes us confident of the validity of eqs. (56,57) in the whole parameter space. In the 
following we will use Landau-Zener formulas to obtain simple analytic predictions for 
the whole range of masses and mixing angles. 

2. As also noted by [19], momentum dependence disappears from eqs. (56,57) (and 
(53,54)). This means that resonance does not produce, as a net result ^, a strong 
spectral distortion of neutrino distributions. As a consequence, the probabilities Paa 
and Pas can be taken, in first approximation, as global multiplicative factors of active 
and sterile neutrino distribution functions. 

So far we have considered neutrino evolution through resonance. Post resonance evolu- 
tion is simple if the mixing angle is sufficiently small. It this case, after resonance, sterile 
neutrinos are not produced. We have then: 

Pss{x,y) _ _ ^Qj^g^ ^ ^ ^^^^ j^gg-j 



Paa{x,y)i 



On the other hand, active neutrino distribution function evolves due annihilation and/or 

scattering, as described by the collision integral Icouiq)- A simple description can be ob- 
tained if active neutrinos are approximately in kinetic equilibrium. Under this assumption, 
the time evolution of the normalization factor is described by the equation: 

^^ = -(7ann)[n„(x)2-l] (60) 

where the symbol (...) indicates averaging over the Maxwell-Boltzmann distribution and 
7ann IS the total I'a'^a annihilation rate. If we assume that the deviation from chemical 
equilibrium produced by resonance is generated abruptly at the time Xres when neutrinos 



^We remark that neutrinos with different momenta go through resonance at different times, sec eq. (46). 
This means that, while we do not expect strong distortions at the end of the process (i.e. after all relevant 
momenta have passed through resonance), a strong distortions can be generated at a given time as a result 
of the fact that only some momenta, at that time, have reached resonance. 



12 



with momentum y = 3 (which we take representative of the whole ensemble) goes through 
resonance, the previous equation can be solved explicitly giving the result: 

where 



P — 1 
A = — exp 

p -\- ^ 



o \7ann/ 1 3 3 



(62) 



and Xggp^ is the "time" at which BBN takes place. 

In figs. 3 and 4 we compare the analytical estimates with the numerical calculations. 
Specifically, in fig. 3 we compare, in the case of v^. — Vg mixing, the prediction of eqs. (57,59) 
(red dotted lines) with the energy density of sterile neutrinos obtained from numerical 
calculations (solid lines) as functions of the mixing angle sin^ 29, for selected values of the 
mass difference 5m?. One can see that eqs. (57,59) correctly describe the dependence of 
sterile neutrino production on the oscillation parameters over a wide range of masses and 
mixing angles. In particular, the analytical calculations are very accurate at small mixings, 
while they tend to overestimate sterile neutrino production at large angles. We remark 
that numerical results shown in fig. 3 have been obtained with the complete code (i.e. 
without using stationary point and kinetic equilibrium approximations) in order to be sure 
not to introduce the errors of approximation into the calculation. In fig. 4 we compare 
the analytical results (dotted lines) for ng with the energy density of Ve. obtained from 
numerical calculations (solid lines). Blue dotted lines are obtained from eqs. (57) and do 
not take into account the post-resonance evolution of Vf>. Red dotted lines include post 
resonance evolution according to eq. (61). As one can see, analytic results give a correct 
qualitative understanding of ng dependence on oscillation parameters. 

4.3 Resonance case. Numerical calculations 

Stationary point approximation (34) is also valid in the resonance case but the region of va- 
lidity is found in a different way because the effective production temperature is determined 
by the resonance condition and not by eq. (24). 

In order to discuss the validity of the stationary point approximation, we have to go 
back to eq. (33). The stationary point approximation follows from the expansion: 



i J^^ dx2- 



of the exponent in eq. (33). Due to the fact that jas/H ^ 1, only a small region Sx = 
{x — xi)/xi ~ H/^ias is important. Thus the expansion can be truncated to first order and 
the stationary point approximation can be easily obtained: 

^""^ (Paa - Pss) (64) 



(Haa -Hss) - H 



In resonance case the situation is more complicated. Far from the resonance, the above 
discussion applies without any modification. At the resonance point, instead, the factor 
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(^aa ~ ^ss) vanishes, with the consequence that higher order terms in the expansion may 
become important. This does not happen when: 



diHa 



n. 



dx 



/ res 



(65) 



where dx = (H/jasJres- In this case, which corresponds drn^ > 10~^ eV^ for — Ug mixing 
and Sni^ > 10~^ eV^ for — Vg mixing, the first order term in the expansion is anyhow 
dominant and eq. (64) continues to be vahd. 

For smaller masses, eq. (64) is not strictly true. Specifically, this equation overestimates 
Pas at resonance. However, this does not introduce significant errors for the evolution 
of the diagonal components, since at the same time, resonance width is underestimated 
by a compensating amount. The most direct way to understand this is to consider the 
equations describing the evolution of the diagonal components which are obtained within 
the stationary point approximation: 



Hxd^{paa - Pss) 



[Ha 



{Paa - Pss) - Icoll{p) 



(66) 



Close to resonance, we can neglect the collision integral and expand the r.h.s obtaining: 

, {jas/HxUinas/Hx)l 



dx{Paa - Pss) 



B{x - Xres)^ + ilas/Hx)'^^ 



{Paa - Pss) 



where 



B = 



\d{Haa-nss)/Hx-\ 


2± —i 


'd''F{x,y) 


dx 


res 


dx"^ 



(67) 



(68) 



The index [res] indicates that the various quantities are evaluated at the resonance point 
and the function F(x,y) is defined in the previous section by eq. (43). If we define ^ = 
B{'yas/Hx)~^{x — Xres), eq. (67) becomes 



diiPa 



(Has/Hx) 



2 

res 



B 



{pa 



Ps 



(69) 



One should note that the factor has disappeared from our formulas. Wc could expect 
then that eq. (66) , derived under assumption of large 7^3 , is valid even in the limit of small 
7as. We can confirm this by integrating eq. (66) assuming a slowly variation of paa — pss- 
We obtain: 



A(/9ao - Pss) - -^-^{Tias/Hx) 



2 

res 



dx^ 



{Paa Pss) I 



(70) 



This result coincides with the expected behavior, eq. (51), and confirm the validity of the 
stationary point approximation in the resonance case even for small mass differences. 

In fig. 5 {Ufj^ — Vs mixing) and in fig. 6 {vg — Ps mixing) we show our numerical results 
for the resonance case in the plane (sin^ 29, 6m'^). Let us remark that, besides the previous 
discussion, we also have checked numerically the displayed results, by comparing them, 
for several values of the parameters, with the results of complete numerical integration of 
kinetic equations. Similarly to fig. 1 and fig. 2 , the first panel shows the effect of oscillations 
on Ug energy density, the second panel shows the energy density of sterile neutrinos, the 
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third panel describes the total neutrino contribution to the energy density and, finally, the 
fourth panel shows the total effects of oscillations on BBN. 

Our constraints on sterile neutrino admixture in the resonant case are stronger than 
previous results of [19], which were obtained in the context of single momentum approxima- 
tion We are not able to understand the source of the observed difference. We are however 
confident in the results of our numerical calculations, especially in the small mixing angles 
region, where, as shown in the previous section, they can be analytically understood with 
good accuracy. 

In fig. 5 and 6 we show with red dotted lines the prediction obtained from eq. (57), 
which correctly reproduce (at small mixing) the numerical results. In Fig. 6 wc also show 
the prediction which can be obtained by using eqs. (56) and (61) which give a qualitative 
understanding of the behavior observed for AA^^^^. 



5 Mixed two active and one sterile neutrinos 

In the case that only two active neutrinos (say, v^. and v^^ are mixed between themselves 
and one sterile neutrino, all algebraic manipulations are relatively easy and the resulting 
equations have a simple form and can be explicitly written down below. We will consider 
large mass difi^erences, so for the off-diagonal elements of the density matrix the stationary 
point approximation may be used. 

Using eq. (7) we can eliminate the e/x-off-diagonal component of the neutrino density 
matrix: 

^ ^e/u(P/U/i Pee) "I" "^esPsii, "^jxsPes (71) 

Substituting this expression into eqs. (8) and their complex conjugate we obtain a closed 

system of equations for four off-diagonal component pse, Ps^j and their complex conjugate 
Pes 

and p^s- They can be explicitly solved in terms of p^^, Peo 

and pss- For example, 

Pse = {AC* - BC) /Det (72) 

where Det = \A\'^ — \B\'^ \s the determinant of the linear system of 4 algebraic equations for 
Psa and Pas and 

A = Ks^^- r..:: '1 2 + — / .X (73) 



B = 2Im 



'Hfj^s'Hes'Hefj, 
Zjj,sZ^^ T~ies 



(74) 



C = Se ^ , Sa + , S* (75) 

^ixs-^jxe ''-es ^fis^fie ''-es 

Zaf^ = 'H^al-'H^f^l+naP (76) 

= {pnn - Pee) + Hes (Pee - Pss) (77) 

Sfj, and psfj, can be obtained from eqs. (71)- (77) by the interchange of e and p. 

Analytical expressions are particularly simple in the case of weak mixing with and 

for sufficiently large mass differences, such that refilling of active states is fast and their 
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number densities approach equilibrium values, Pee = P/i^^ ~ /eg- The analytical limits (38) 
and (39) are obtained under similar assumptions. In this case the differential equation that 
governs evolution of Pss takes the form: 

dn t r \ Tes (^es^/i/i ^e/i^/xs) "I" Tm* (^/is^ee ^e/i^es) fno\ 
xPss = ^ [Jeq - Pss) 72 \'°) 

Hx [TCee'HfMfj, — Tien ) 

This result is valid in no-resonance case when Tiee^ufi ~'^eii 7^ 0- 1^ MSW-rcsonance exists 
then one should take into account imaginary parts proportional to coherence breaking terms 
7. In numerical calculations such terms have been included. 

Simple expression (78) has been derived under assumption of weak mixing between 
active and sterile neutrinos. Relaxing this approximation we find that resonance condition 
is: 

A^eAe.A,^ - n%A^s - nlAse - H^Ae^, = (79) 

where Aq,^ = TCaa - Upfs. 

Eq. (78) justifies neglecting of non-diagonal terms in the effective potential. Indeed 
these non-diagonal terms are proportional to the non-diagonal components of the active- 
active part of neutrino density matrix. As we see from eq. (78) they should be compared 
to the diagonal ones. In the situation close to thermal equilibrium p^b = if a 7^ 6, while 

Paa = /eg ^ 1- 



6 Realistic case of three active and one sterile neutrinos 

Recent data from SNO [1], KamLAND [2], and earlier data on solar [3] and atmospheric [4] 
neutrinos have provided strong indications in favor of neutrino oscillations. It is practically 
established now that all active neutrinos are mixed with parameters given by the Large 
Mixing Angle (LMA) solution to solar neutrino problem [11], see eq. (1) and by atmospheric 
neutrino data [4], eq. (2). Existence of fast transitions between fg, f^, and Vj- may change 
BBN bound on mixing with sterile neutrinos especially for small values of mass difference. 
This means that one has to take into account mutual active neutrino mixing in the discussion 
of BBN bounds. Clearly the situation become much more complicated both conceptually 
and numerically. In the following, to keep the discussion as simple as possible, we will make 
the assumption that active-sterile neutrino mixing is small. 

In the case of 4 mixed neutrinos the algebra becomes significantly more complicated. 
This is why we will not present explicit expressions for off-diagonal components of density 
matrix through the diagonal ones. All this work has been done by computer and the results 
have been inserted into differential equations for diagonal components. The latter have 
been solved numerically. To check the accuracy of the procedure we solved numerically 
the complete system of kinetic equations (5)- (8) for certain fixed values of the parameters 
without any simplification and found that agreement with the approximate numerical results 
is very good. We did not do that for all interesting values of the parameter because it would 
take too much computer time, especially in the resonance case. 

It is clear that, depending on masses and mixing, one can expect in general case one or 
more resonances among active and sterile neutrinos. If 1^4 is the heaviest mass eigenstate 
there may only be resonances in the transformations among active neutrinos (they are 
quite efficient anyhow because of the strong mixing). As is becoming lighter and lighter, 
there appear first one resonance in Va-Vg transformation, then two and ultimately three 
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resonances. The resonance condition is given by the vanishing of the determinant of the 
6 X 6-matrix: 



Dq = det 
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(80) 



In the limit of weak mixing with i^s the determinant of this matrix is reduced to the product 
of two determinants: 
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(81) 



The first term determines resonances in the active neutrino sector, while the second one 
determines resonances between Vg and three active neutrinos. Of course resonances manifest 
themselves in the process of numerical calculations and we have checked that they are indeed 
in the right positions determined by the above equations. 

6.1 Simplifying the problem 

If we assume that active-sterile neutrino mixing is small, we can write the neutrino mixing 
matrix as: 



U 



Tract 
rract 



rract 
rract 

n'i 

Tract 
'-'t2 

£2 



Tract 

'-'ns 

Tract 
£3 



m \ 
1 ) 



(82) 



where = U23 ■ f/ia • U12 is the 3x3 active neutrino mixing matrix and: 



(83) 
(84) 
(85) 



•'atmo- 



The previous relations are correct to first order in £1,62,63 (or equivalently /71, 772, ^73). 

We assume for simplicity that ^13 = 0, while we take Ou = ^soi and ^23 
Concerning mass differences, we assume 6m^2i = "t-I ~ "^1 = Sm^soij and 5m^32 = 
1712 = Sm^atmo- Solar and atmospheric oscillation parameters are given by eq. (1) and (2) 
respectively. 

If ^13 = 0, one can take advantage of the symmetry between and Vt to simplify the 
problem. The early universe, indeed, does not distinguish between and z^t- at T < m^. 
We can thus make a rotation in the flavor basis: 



I'. 



23 



/ ye \ 



(86) 
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In this basis the mixing matrix becomes: 

I cos 6112 
- sin 6*12 



where: 



U' = 



V 





ei 



sin 6*12 

cos 9\2 



£2 





1 

£3 



i2 



V's 



(87) 



= m = -cos6'i2 



COS 023 ■ 112 — sin 023 ' % 
sin 023 • ??2 + COS 023 " ^ 



ei — sin 012 • £2 

= sin 012 • ei — cos 0i2 • €2 



= £3 



(88) 
(89) 
(90) 



The problem which was written in terms of the parameters (r/i, 772, 773) can be easier treated 
in terms of the parameters (^/i • ??2' "^/s ~ ^s)- using these parameters, it can be easily 
separated into two smaller sub-problems. In the following sections, we will discuss separately 
the constraint on 773 and (??i,r72)- In the final section we will combine our results. 



6.2 Bounds on 7/3 

The bounds on 773 = es can be immediately obtained from the results of section 4, since this 
is a pure "1 active + 1 sterile neutrino" problem. Clearly the result depends on the value 
of the mass difference 

5m\^ = m\ — m\ . (91) 

When 6m1^ > active-sterile transitions are non-resonant. As a consequence, eq. (39) (and 
subsequent discussion) is valid. In terms of rj'^ we have thus: 

(5m243/eV2) (2?7^)4 = 1.74 • lO'^ ln\l - AN^) (92) 

When (^777.43 < 0, sterile neutrino production occurs though resonant transitions. This 
means that eq. (57) (and subsequent discussion) is approximatively valid. With very simple 
algebra, assuming that the number of extra neutrinos AA^,^ is approximatively given by the 
sterile neutrino energy density, we can obtain the following relation: 

{\Smh3\/eV'^) i^v's)^ = 1-9 • 10"'' 111^(1 - AN,,) (93) 

We remark that the two equations above take into account only sterile neutrino production, 
while they do not consider the effects related to depletion of other neutrino flavors. As 
a consequence, they are quite accurate for large mass differences, while they generally 
underestimate the effect on BBN at smaller masses (since they neglect effects of depletion 
on n — p transformations) . More precise bounds can be read from fig. 1 and 5 with the 
identification Sm^ — Sm^^ and sin^(20) — (2773)^. 



6.3 Bounds on 'r][ and 773 

The problem described by the parameters (771, 772) is much more complicated since it involves 
directly two active and one sterile neutrino, with mixing described by: 
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The mixing angle O12 and the mass difference 617121 fixed by KamLAND and by solar 
neutrino experiments, see eqs. (1). 

It is very difficult to obtain analytic estimates in the general case. However, one can 
simply understand what happen in the limiting cases. Let us consider the situation 1114 3> 
1712 > 1^1 (or rui <C mi < m2). In this case, we can approximatively neglect the mass 
difference (5m^2i- The angle O12 can thus be rotated away, resulting in: 
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(95) 



By repeating the discussion of the previous section, one obtains then: 

((5m24i/eV2) {2r][f = 3.16 • 10"^ ln2(l - AA^^) (96) 
(5m242/eV2) (2r?^)^ = 1.74 • 10"^ ln2(l - AA^^) (97) 

where we assumed positive mass differences, i.e. 5m^4i ~ 6171^42 S> Sm^2i > 0- When mass 
differences are negative, i.e. 6171^41 ~ Sm?42 ^ —Sm^2i < 0, one expect instead 

{\Sm'^ 4i\/eV^) {2 T][)^ = 5.2 ■ 10"^° ln2(l - AA^) (98) 
{\Sm\2\/eV'^) {27)2)'^ = 1.9 • IQ-^ ln2(l - AAT^) (99) 

The other limiting situation that we can understand analytically is 777,2 3> 7714 ~ mi. 
In this case we can approximately neglect the mass difference 6171^41. The mixing in the 
(1 — 4)-sector can thus be re-absorbed, leading to: 
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where: 

(f = ri[- tan ^12 + ri2 = "^2/ cos(6'i2) (101) 

We can see that the mixing between active and sterile neutrinos is described by a parameter 
if. If 7^ 0, a resonant transitions occur among Vs and v'^. This mechanism should dominate 
the production of sterile neutrinos. As a consequence, we expect the following relation to 
be valid: 

(|(5m242|/eV2) (2(^-cos(^i2))^ = 1.9 • 10"^ ln2(l - AAT^) (102) 

where (p is given as a function of rj'i and rf2 by eq. (101) ^. 

In order to check the validity of the previous conclusions and to have a better under- 
standing of the BBN bounds on sterile neutrinos, we solved numerically kinetic equations, 
varying, one at a time, the parameter r][ and 772- As a first step, we considered the situation 

®We note that the parameters rj'i and 772 can be chosen in principle in such a way that ip = 0. When 
if = 0, sterile neutrinos are not mixed with active ones and we should expect that sterile neutrinos are not 
produced. In reality, this result simply indicates that, in the limit considered (m2 3> m4 ^ mi), no eflFects 
driven by the dominant mass difference, 5m^2i are expected. Even in this case, however, sterile neutrinos 
can in principle be produced, due to oscillation driven by sub-dominant mass differences 6m^ 41 ■ 
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when the mixing matrix U can be written as U = C/24 ' U12, with 8w?2\ and d\2 fixed at 
solar values, see eq. (1). It is evident that, in the limit of small mixing, one has ^24 — ?72- 
The results obtained are shown, as functions of 5m^42 and sin^(^24)) in figs. 7 (for the case 
5m^42 > 0) and in figs. 8 (for the case (5m^42 < 0). Red dotted lines correspond to analytic 
estimates which are obtained from eq. (97) in the case of positive (5m^42 and from eq. (99) 
in the case of negative 5m^42. One can see that there is a good agreement with numerical 
calculations for large values of mass difference. For small values of |(5m^42| agreement is less 
satisfactory (especially for 5m^42 < 0) both because the analytic estimates, derived for the 
limiting behaviours, are less appropriate and because depopulation of electron neutrinos 
becomes relevant. In particular, by comparing figs. 7 and 8 with figs. 1 and 5 (which show 
the results obtained when we neglect the mixing J7i2), one can see that generally the effect 
of mixing between active neutrinos is to increase v^. depletion, with the consequence of a 
stronger bound on the active-sterile mixing. 

Similar conclusions can be obtained from figs. 9 and 10. In these figures we show 
numerical results for the situation in which the mixing matrix U can be written as C/ = 
t^i4 • V\2 . In the limit of small mixing angle one has ~ r/i . Analytic estimates (red dotted 
lines) are obtained from cq. (96) for 5m? 41 > drn^^i, from cq. (102) for < Sm^^i < 5m? 21 
and from eq. (98) for 5m^ 41 < 0. One sees that there is a good agreement between numerical 
calculations and analytical estimates. In all the cases, analytic estimates underestimate the 
BBN bounds on active sterile mixing. 

6.4 Combining results 

As a summary, by considering variations of the parameter 773 = 63, we have obtained: 

{6m^ 43 / ey^) i'^Va)'^ = 1.74 ■ 10"^ ln2(l - AiV^) for 7714 > mg (103) 
(|(5m243|/eV2) {2ri's)'^ = 1.9 • lO"'^ ln'^{l - AN^) for m4 < mg (104) 

These relations are approximate and generally underestimate the effect on BBN at small 
masses. More precise bounds can be read from fig. 5 with the identification Sm"^ — ^ ^1^43 
and sin2(2^) (2r?^)2. 

By considering variations, one at a time, of rj'i and r]2 we have obtained: 

3.16 • 10"^ln2(l - AiV^) 
1.74- 10"'Mn2(l - AN^) 

for 777-4 > 7772 > 777l (105) 



(W4l/cV2) (2 7?0' = 

(W42/eV2) (27?^)^ = 



(15^^42 l/eV^) [2 (r/isin(0i2)+r/^ cos 012)]^ = 1.9 • 10"^ ln2(l - AAT^) 

for 77l2 > 7714 > 7771 (106) 

(1(5777241 l/eV^) (2?7i)^ = 5.2 • 10-i°ln2(l - AN^) 
(1^777242 |/eV2) (27/^)^ = 1.9 • 10-9 ln2(l-AiV^) 

for 77l2 > 7771 > 7774 (107) 

As above, these constraints are approximate and more precise bounds can be read from figs. 
7,8,9,10. 

In general scenario all parameters are possibly different from zero at the same time and 
the effects described above are all observed togheter. We do not expect, however, strong 
compensations to occur because the described effects have generally different magnitude 
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and because, in terms of the parameter , ??2' ^3)5 problem is well separated in nearly 
all range of possible masses (except for m2 > m4 > mi). This means that the bounds on 
(77^, r/2, ^^3) can be translated into constraints on (?7i, 7721 %) and (ei,e2,e3). We remind that 
V = U23- v' and e = -(C/12)"^ • rj'. 

7 Non-vanishing lepton asymmetry 

Wc assumed above that lepton asymmetry of cosmological plasma is cither very small or 
has the natural value, rj ~ 10"^ of the same order of magnitude as cosmological baryon 
asymmetry. The charge asymmetric contribution to neutrino effective potential in plasma 
(eq. (10)) is 

y(as2/m) ^ -^ -^ . 10-20 (^(«)/10-9)/x^ (108) 

Thus one can see that for such r] the Hamiltonian (9) is either dominated by the mass dif- 
ference term (at low temperatures T = l/x) or by the non-local one (at high temperatures). 
However the magnitude of lepton asymmetry of neutrinos can be much larger than 10"^ 
either due to resonance rise of the asymmetry because of charge asymmetric transformation 
of Ua into Ug and v'a into i>s [22, 33] or due to primordially generated large lepton asymmetry. 
There are plenty cosmological scenarios leading simultaneously to small baryon asymmetry 
and to much larger, even of order unity, lepton asymmetry [34]. In both cases neglected up 
to this point charge asymmetric contribution into Vaa (10) could become dominant. 

The best bounds on the values of the lepton asymmetry can be found from BBN, see 
review [20]. These bounds have been considerably improved due to strong mixing between 
active neutrinos [27, 35]. It allows to put the upper bound on the common chemical potential 
of all neutrino species 

l^al = \l^a/T\ < 0.07 (109) 

This limit is loose enough to allow the charge asymmetric part of the potential to play 
a significant role in neutrino oscillations in the early universe. An analysis of BBN with 
non-zero chemical potentials and possible active-sterile neutrino mixing has been performed 
recently in ref. [21]. 

Let us first consider naturally small initial value oirj < 10^^ and later more questionable 
larger values of primordial lepton asymmetry. For small asymmetry and positive mass 
difference between Ug and any active neutrino MSW-resonance would not be efficient (see 
below eq. (122). Hence lepton asymmetry would remain always small and (108) 
may be neglected. Thus the bounds obtained above in non-resonance case are justified. 
This may be not so if mass difference is negative and the resonance transition between Ug 
and Ua could take place. As a result a large lepton asymmetry, close to or even a little larger 
than rj = 0.1, can be generated in the active neutrino sector [22, 33]. Of course the total 
leptonic charge of active plus sterile neutrinos is conserved and a lepton asymmetry of active 
neutrinos is compensated by opposite sign of lepton asymmetry of Ug. Neglecting y(««2/"») 
in this case would be unreasonable. However, as we see below, a large lepton asymmetry 
can be generated only for very small values of mixing angle between Vg and i>a- 

sm^2^|5m2|<| 3^^Q_, J^^^ (110) 

Here mass difference is measured in eV^. Note that the bound for is stronger than for f/i,r) 
in contrast to the bounds (38), (39). For the values of parameters outside this region the 
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approach used above (with neglected) is justified. As we have seen in the previous 

sections, the allowed interval of {O—Sw?) found in the resonance case but without ^ is 

inside the parameter space restricted by (110). Hence we can conclude that the area outside 
of (110) is definitely excluded and that the real bound can be stronger. It is more difficult 
to make conclusion about allowed range of the mixing parameters inside the region (110) 
because an excessive number of neutrino species may be compensated by electronic charge 
asymmetry. However, one should take into account that the generated lepton asymmetry 
could be strongly inhomogeneous [36] (see also [20, 37]). If the mixing angle is larger than 
10"^, so that the lepton asymmetry in the minimum (see below discussion after eq. (117)) 
drops more than by 5 orders of magnitude, then the universe would consist of chaotic 
domains with cqTial positive and negative values of ry [20] and the impact of the electronic 
lepton asymmetry on BBN would be less pronounced after averaging over such domains. 
On the other hand, as argued in ref. [37], production of sterile neutrinos on very steep 
domain walls might be resonance enhanced and it may lead to unacceptable large AN^,. 
The problem is not yet settled down and deserves more investigation. 

Anyhow, even the limit (110), which allows the resonance generation of lepton asym- 
metry is quite strong, much stronger than the non-resonant one and, for |(5m^| > 10~^, it 
is also stronger than the limits obtained in direct measurements. 

Resonance rise of lepton asymmetry does not take place in the case of strong (vg-i^a)- 
mixing by the following evident reason. If mixing is large both active and sterile neutrinos 
quickly reach thermal equilibrium when charge asymmetry is vanishing, paa = Paa = Pss = 
Pss = feq- To find the magnitude of mixing angle when fast equilibration takes place we will 
use analytical calculations of ref. [33]. The evolution of the charge asymmetry is governed 
by the following equation derived in the limit of large values of parameters I: 



dZa{q) lQ^^Kl{sm2ef 



dq 



167r2 



POO 

/ dyfegiy) 

Jo 



/ dqiT.{qi,y) exp 

Qin 



qqi 



sm 



y qqi 



sm 



baKa I dq 

Jqi 



Za{q2) 
' 4/3 
12 



(111) 



where ( = Ka{q — qi) is effectively integration variable and the constants depending upon 
the active neutrino flavor are given by: 



Ke = 5.63 -10^(1^1 cos 2^)^/^ K^,. 



6, = 3.3- 10-3(|(5m2| cos2^)-^/3^ 
«i 7.4•10~^ 



-'e 



b 



fJ.,T 



-- 2.97- 10^(1 (5m2 1 cos 2^) 
7.8- 10-3(|(5m2| cos2e) 
i 5.2 • 10"^ 



1/2 
-1/3 



(112) 
(113) 
(114) 



The quantity Z is related to charge asymmetry of active neutrinos as 

Za = 1.67 • lO^a 



(115) 



and the new "time" variable q is chosen in such a way that in the limit of small charge 
asymmetry the resonance always takes place aX q = y, i.e. q = PaX^ with 



1 /2 

/3e = 6.63 • 10^ (|(5mVos26l) , /3^,r = 1.257 • 10*^ (|<5mVos 26* 
The function S(g, y) is given by 

^ = Paa + paa ~ Pss ~ Pss 



1/2 



(116) 



(117) 
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Equation (111) was solved in ref. [33] under assumption that S varies very slowly and 
in an essential interval of q it is approximately equal to its initial value, Sj„ = 2. This 

approximation is valid for sufficiently small mixing angle 6 and we will determine below 
when it is true. For small q the charge asymmetry Z exponentially decreases and only for 
q > 1.278 it starts to rise. For sin 20 = 10~^ and drn^ = 1 eV^ the asymmetry drops by 3 
orders of magnitude, for larger 6 the decrease would be much stronger. 

In the limit of small charge asymmetry S satisfies the equation (in this equation and 
below we will omit sub-indices a): 



dq \ y 



dqi^{qi,y) exp 















cos 






- qqi. 






qqiJi 



(118) 



For large K the integral in the r.h.s. "sits" on the upper limit and can be easily taken. Thus 
we come to the equation: 

dT, Ke II cr siu^ 20 „ , , 

dq [q^ - y^) + e^y* 

The solution of this equation is straightforward. It can be substituted into eq. (Ill) and 
we find that asymmetry does not rise if E quickly decreases. According to this solution, 
this happens if the mixing is bounded from below as given by eq. (110). 

Let us now consider the case of a large primordial lepton asymmetry and check if a 
stronger mixing between active and sterile neutrinos is allowed. If so, then an experimental 
discovery of mixing between violating the bounds derived above would indicate that 
the cosmological lepton asymmetry is much larger then the baryonic one or that neutrinos 
possess a new interaction, as e.g. coupling to Majorons, which might strongly change the 
oscillation behavior in cosmological plasma [38]. 

If the contribution of the lepton asymmetry into neutrino effective potential (10) is 
non-negligible, then resonance transition would exist for any sign of mass difference. In 
both cases, there are two resonances but for SmP' > these two resonances are either in 
neutrino or antineutrino channels (depending on the sign of rj), while for 6m'^ < there is 
one resonance in neutrino channel and one resonance in antineutrino channel. We assume 
that the mixing is sufficiently weak, so only the resonance transitions are of importance. 
The resonance condition looks as following: 

5 • W-^^C2 6m^ - + 10-20 «;„ ^ ± 1.1 • lO'^o j^g J- = (120) 
y x"^ 

where C2 = cos 20, Kg = 1-137, /t/i^r = 0.317, r^g = lO^r/, the signs ± refer to antineutrino 
(-I-) and neutrino (— ) respectively. Solving this equation we find the resonance values of 
neutrino momentum: 

yi,2 = ±(0.55/K)r/9 (l ± (121) 

where A = 1.65 • l{f'C2 K.6m? /rj^. 

If dm? > 0, i.e. Us is heavier than active neutrinos, then the resonances occur only for 
Xx^ < 1. The oscillations could have an impact on BBN if they are efficient in the interval 
X = 0.1 — 2. Hence the charge asymmetry should be sufficiently large: 

%>1.3-10^c^/'acV2(<5^2)1/2 (122) 
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We can roughly estimate resonance contribution into production of sterile neutrinos 
integrating eq. (36) over y with pss < Paa = feq- 




2A5 ■ 10''' K-'sic2{Sm^y X 



2x2 8 yi exp(-yi) + yi exp(-yi) 



yi - y2 



(123) 



where S2 = sin 20. Equation (36) gives a good description of the production rate also in 
the resonance case if the mass difference is not too small. As we have mentioned above, both 
resonances are effective either in (fa-z/s)-transition (for > 0) or in (t'a-^'5)-transition (for 
T] < 0). Hence resonance transitions diminish the absolute value of the original asymmetry. 
Indeed, if e.g. the resonances operate in Va ^ channel the transition of Va into Vg would be 
more efficient than the charge conjugated transition of Va into Vg and the number difference 
between Va and Ua would decrease. If the resonance took place above T = 3 — 5 MeV (or 
at X below x = 0.33 — 0.2), then the distribution of Ua and Ua would have equilibrium form 
with chemical potentials satisfying = —^a and the standard BBN code can be applied 
to calculations of light element production. For smaller resonance temperatures chemical 
equilibrium would not be restored and thus 7^ —^a- The code should be modified to 
include this effect. 

Let us find now what is the largest possible value of the lepton charge asymmetry r] 
allowed by the limit (109): 



3 vrV6 
11 1.8 ^ 



0.25^ 



(124) 



where the first factor 3/11 comes from the present day neutrino-to-photon number ratio and 
the second one comes from the ratio of the equilibrium number densities {n, — np) / {uu +np). 
The result is valid in the limit of small ^. Since individual chemical potentials before 
equilibration enforced by v^ — v^ — Vr oscillations may be at most three times larger than the 
limiting value |^| = 0.07 we conclude that the maximum magnitude of the primordial lepton 
asymmetry is r/ = 0.05. Another formally open possibility is large values of 2r]y^ + rji,^ + rjn^ 
and/or of "Irji,^ + r]y^ + rji,^ (and fi ^ t), contributing into effective potential of different 
neutrino species but a strong compensation in the sum + + down to 0.2. 

If we take rj = 0.05 then the contribution of the first resonance is strongly suppressed 
and only the second resonance corresponding to smaller y is effective. The produced Vg 
would contribute to the additional number of neutrino species as: 



AN, = 0.8 • 10^sl{6m^ 



dxx^^ exp 



-0.9(5m^a;^ 



(125) 



(e) _ 



0.3 in 



where it was assumed that C2 = 1. As an upper limit of integration we will take x. 
accordance with the estimates presented after eq. (17). If we mixing would be predominantly 



with or V. 



one should seemingly take Xmax = 0.2. However this is not so because 
refilling of active neutrino species would take place through e+ v^. with subsequent 

fast transformation of into and through oscillations and hence one should take 

™(/^'^) _ _ n 

Xmax — ■i^max — u.o. 

If the mass difference is smaller than 5m? max ~ 10^ eV^ (notice that cosmological limits 
on masses of active neutrinos [20] are not applicable to Vg because their number density 
may be much smaller than the equilibrium one) then exp[— 0.9(5m^a;^] 1 and we obtain 



AN, = 7 • lO^s^ ( W 



11 



lO^s^ {5rr? 



(126) 
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If S2 ~ 0.1 then the oscillations could diminish primordial lepton asymmetry of active 
neutrinos and produce noticeable amount of Vg- 

However a more important for BBN effect comes from larger x. Above x = 0.3 refilling of 
neutrino states may be neglected and the total number/energy density of + f ^ + + 
would be constant. Correspondingly production of sterile neutrino states would lead to 
an equal decrease of active ones. Due to strong mixing between the later the deficit of 
active neutrinos would be equally distributed between all three active species. As we have 
discussed above, BBN is very sensitive to a deficit of V(>- We estimate the Vs production 
using eq. (125) with Xmax = 1-4, i.e. inverse temperature of neutron-proton freezing. For 
dm? = 1 eV^ we obtain AN^^ = 3 x lO^s^ and for small drn^, such that 6m'^ x^ < 1, 
A.N,y^ = 3 X 10^ S2{Sm^)^ . The deficit of v^, (or 1?^, depending on the sign of initial ij) would 
be equal to one third of that. However there is one more twist, namely there was some 
primordial lepton asymmetry of active neutrinos and oscillations between and z^^ would 
diminish the absolute magnitude of the latter. As a result we could arrive to (almost) equal 
number densities of active neutrinos and antincutrinos which would be suppressed with 

respect to the standard one by the factor exp[— | '^a 1/3] where Ca^"*^ are primordial 
values of chemical potentials of active neutrinos. According to eq. (28), a deficit of Ve 
and equal deficit of v'e is equivalent to addition of AN^ = 12.3|^|. In this way we obtain 
essentially the same bound as eq. (109). However, it is not so in a general case. Invoking 
an additional parameter, the value of primordial lepton asymmetry, evidently makes the 
bounds on mixing between and Ua less restrictive as well as the bounds on the value of 
lepton asymmetry obtained from BBN in absence of mixing with Vg- 

Now let us consider negative mass difference, 5m? < 0, when resonance is possible even 
with vanishingly small charge asymmetry. Since active neutrinos are assumed to be heavier 
than Ug the mass difference is bounded by the recent WMAP data [40] by drn^ < 0.06 eV^. 
In this case, of negative Sm^, there are again two resonances but now one is in neutrino 
channel, while the other is in antineutrino channel, as can be seen from eq. (121) with A < 0. 
If the initial value of \r]\ is sufficiently large, then \r)\ would rise as a result of oscillations 
between !/„ and Ug in contrast to positive mass difference, dm^ > 0, when oscillations lead 
to a decrease of |?7|. 

The rate of transformation of Ua into Ug and of the charged conjugated process can be 
found from equation which is similar to eq. (123) where only one resonance is effective. 
For a positive lepton asymmetry, ?? > 0, the resonance in neutrino channel takes place at 
y = yi (with positive sign of the square root) and the resonance in antineutrino channel 
takes place at y = y2- For a large and, say, positive r] the resonance in neutrino channel is 
too far, yi ^ 1, in the essential range of x and its contribution can be neglected - this is 
similar to the discussed above case of positive 6m^. The resonance in antineutrino channel 
takes place at y2 ~ (4.5 • 10^/779)2;^ and could efficiently transform i>a into Ug. It would lead 
to an increase of rj. For a more accurate calculation of the transformation rate one should 
include non- vanishing chemical potential into number density of active neutrinos in eq. (36), 
Paa = feq = exp(— y + ^) but since we usually consider small ^ <C 1, this correction is not 
of much importance. The oscillations create an increase of the absolute value of the lepton 
asymmetry and hence the BBN bounds on the initial magnitude of the latter should be 
generally stronger than in no-oscillation case. However there are several effects "working" 
in the opposite direction which may compensate each other and so the bounds on r] for 
concrete values of the oscillation parameters in each case should be derived separately. It 
can be approximately done using results presented in this section. 



25 



8 Conclusion 



Wc have shown above that the effects of strong mixing between all three active neutrinos 
noticeably change the BBN bounds on the possible admixture of sterile neutrinos previously 
obtained under assumption of no active-active mixing. Our results are summarized in sec. 6 
with simple analytical expressions. They essentially depend upon the sign of mass difference 
between different neutrino species, because the latter determines the number and positions 
of the resonances. More detailed and more accurate bounds can be read-off from figures 
1,5,7,8,9,10, panel 4 in all the cases. The results depend upon the limits on the number of 
effective neutrino flavors allowed by BBN. One can see that even with very modest limits 
there is not much freedom permitted to admixture of sterile neutrinos to the known active 
ones. 

If lepton asymmetry of the universe is non-negligible the conclusions are less definite 
because of the additional parameter, the value of the asymmetry. However, if asymmetry is 
generated by the oscillations themselves, the mixing with Ug should be very small, below the 
bounds derived for the case of vanishing asymmetry. For the case of primordial asymmetry 
the bounds are weaker and with some conspiracy they could be very weak, allowing rather 
strong mixing between sterile and active neutrinos. More detailed investigation is desirable 
in this case. 
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^f^-i^s mixing - Non resonance case 
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Figure 1: Numerical results for i/^ — Vs mixing. Non resonance case. First panel: Iso- 
countour lines for fe energy density. Each line corresponds to fixed value of z^e energy density. 
The indicated values are normalized to equilibrium value. Second panel: Iso-contour lines 
for Vg energy density. Third panel: Iso-contour lines for total neutrino contribution to the 
energy density. Each line corresponds to a fixed value of the number of extra neutrino flavors 
AA'^^. Fourth panel: Total effect of neutrino oscillation on BBN. Each line corresponds to 
a fixed value of the effective number of neutrinos AA'^^^^ species calculated according to 
eq. (28). Experimental bounds on AA*"^^^ exclude the region above the corresponding 
curve. Red dotted lines correspond to analytical estimates obtained from eq. (39) and 
following discussion. The quantities presented in this and in the following figures are taken 
at the BBN "time", xbbn = 1-4. 
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Figure 2: Numerical results for Ue — i^s mixing, non resonance case. See Fig.l for detailed 
explanation of the various panels and for a definition of the various lines. Red dotted lines 
correspond to analytical estimates obtained from eq.(38) and following discussion. 
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Figure 3: Comparison between numerical results (solid lines) and analytic estimates (red 
dotted lines) for Vf, — mixing, resonance case. The various lines show the energy density 
of sterile neutrinos (at the time of BBN) as a function of the mixing angle sin^(2^), for 
selected values of the mass difference brr? . Analytic estimates are obtained from eq. (57) 
and following discussion. 
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Figure 4: Comparison between numerical results (solid lines) and analytic estimates (dotted 
lines) for V(,—Vs mixing, resonance case. The various lines show the energy density of electron 
neutrinos (at the time of BBN) as a function of the mixing angle sin^(26), for selected values 
of the mass difference 5m?. Blue dotted lines, which are obtained from eqs. (56), do not take 
into account Vf, post-resonance evolution. Red dotted lines take it into account following 
eqs. (61, 62). 
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z^^-i^s mixing - Resononce case 
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Figure 5: Numerical results for — Vg mixing, resonance case. See Fig.l for detailed 
explanation of the various panels and for a definition of the various lines. Red dotted lines 
correspond to analytical estimates obtained from eq.(57) and following discussions. 
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z^e-i^s mixing - Resonance case 
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Figure 6: Numerical results for z/g — Vs mixing, resonance case. See Fig.l for detailed 
explanation of the various panels and for a definition of the various lines. Red dotted lines 
correspond to analytical estimates obtained from eq.(57) and following discussions. 
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Figure 7: Numerical results for the case in which the neutrino mixing matrix can be written 
as U = C/24 • U12, with 6m?2i and 612 fixed according to eq. (1). We consider positive values 
for the mass difference 5m? ^^2 = fn\ — m^- See Fig.l for a detailed explanations of the 
various panels and for a definition of the various lines. Red dotted lines correspond to 
analytic estimates obtained from eq.(97). Note that in the small mixing angle limits, the 
angle 6*24 coincides with the parameter r/2 defined in eq.(89) 
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Figure 8: Numerical results for the case in which the neutrino mixing matrix can be written 
as U = U24 • C/i2) with 5m^2i and 612 fixed according to eq. (1). We consider negative values 
for the mass difference (5m^42 = ml — m^. See Fig.l for detailed explanation of the various 
panels and for a definition of the various lines. Red dotted lines correspond to analytic 
estimates obtained from eq.(99). Note that in the small mixing angle limits, the angle 624 
coincides with the parameter 772 defined in eq.(89) 
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Figure 9: Numerical results for the case in which the neutrino mixing matrix can be written 
as U = Ui4 • Ui2, with 6m?2i and 9i2 fixed according to eq. (1). We consider positive values 
for the mass difference 6m? /^i = — See Fig.l for detailed explanation of the various 
panels and for a definition of the various lines. Red dotted lines correspond to analytic 
estimates obtained from eqs.(96,102). Note that in the small mixing angle limits, the angle 
014 coincides with the parameter 77^ defined in eq.(88) 
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Figure 10: Numerical results for the case in which the neutrino mixing matrix can be written 
as U = Ui4 ■Ui2, with and 612 fixed according to eq. (1). We consider negative values 

for the mass difference (5m^4i = ml — See Fig.l for detailed explanation of the various 
panels and for a definition of the various lines. Red dotted lines correspond to analytic 
estimates obtained from eq.(98). Note that in the small mixing angle limits, the angle 614 
coincides with the parameter -q'^ defined in eq.(88) 
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